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Full-field x-ray microscopy using x-ray objectives has become a mainstay of the bi¬ 
ological and materials sciences. However, the inefficiency of existing objectives at 
x-ray energies above 15 keV has limited the technique to weakly absorbing or two- 
dimensional (2D) samples. Here, we show that significant gains in numerical aperture 
and spatial resolution may be possible at hard x-ray energies by using silicon-based 
optics comprising ‘interdigitated’ refractive silicon lenslets that alternate their fo¬ 
cus between the horizontal and vertical directions. By capitalizing on the nano¬ 
manufacturing processes available to silicon, we show that it is possible to overcome 
the inherent inefficiencies of silicon-based optics and interdigitated geometries. As 
a proof-of-concept of Si-based interdigitated objectives, we demonstrate a prototype 
interdigitated lens with a resolution of ~ 255 nm at 17 keV. 
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X-ray microscopy (XRM) is an established family of techniques for imaging embedded 
and structurally complex specimens with sub-pm resolution. The ability to look inside 
dense matter has provided crucial insight into phenomena such as ferromagnetic domains^, 
nano-scale strain^ and compositional inhomogeneity^. The techniques can be broadly cat¬ 
egorized into scanning-'^, projection-^ and objective-hased^ approaches, which offer different 
compromises of spatial resolution, acquisition time and sensitivity. Full-held imaging with 
an objective is particularly relevant to materials and geological sciences, as its efficiency and 
modalities (e.g. dark-held^) enable real-time imaging of complex processes®. Performing 
full-held XRM at hard x-ray energies (>15 keV) would then open a new door to dynamic, 
three-dimensional (3D) multi-scale studies of denser and more complex samples. However, 
current x-ray objectives tend to be aberrated or inefficient in the hard x-ray regime, limiting 
spatial and temporal resolution. 

In this letter, we show that Si-based 2D compound refractive lenses (CRLs) are a vi¬ 
able approach for improving the numerical aperture (NA) and efficiency of full-held XRM 
objectives at hard energies. Specihcally, we show theoretically that sufficiently miniatur¬ 
ized Si CRLs can outperform the current state-of-the-art at hard energies, and validate this 
prediction based on the performance of a prototype objective at 17 keV. 

Various x-ray imaging optics have been proposed for use at hard energies: rehective 
optics (multilayer mirrors®) are efficient but expensive and delicate; dihractive optics (Fresnel 
zone plates^®) have large NAs but are inefficient, while refractive optics (compound lenses^^ 
or prisms^^) can be efficient but are prone to aberration and small NAs. The approach 
demonstrated here can overcome some of the current limitations of refractive optics at hard 
energies by utilizing a miniaturized, 2D interdigitated conhguration of planar silicon lenses 
(Fig. 1). 

The ideal lenslet geometry and conhguration of a CRL-based XRM can be determined by 
optimizing the NA. From ref. 6, the NA of an XRM using an ideal CRL of focal length /, 
comprising N identical lenslets with apex radius of curvature R, linear attenuation coefficient 
jj, and refractive decrement 5 follows the relationship; 

Eq. 1 implies that for a given magnihcation AT, the NA is greatest in CRLs with a 
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FIG. 1. SEM image of Si interdigitated lenslets within the prototype objective lens tested in this 
work. The small apex radii of curvature accessable to Si-based lenslets enables shorter focal lengths 
and higher NAs than existing 2D objectives 

short focal length and lenslets of a low-Z material. Many CRLs therefore utilize lenslets 
produced by indenting the parabolic lens prohle into polycrystalline Be or Al®. While these 
materials have favourable values of 6 and p in the hard x-ray regime, the indenting process 
is expensive and limited to large R due to material (grain structure, porosity, plasticity) 
and processing (tool shape, concentricity) issues. At hard energies, this necessitates the use 
of many lenslets at signihcant cost to achieve a short / and thus large NA. Polymer CRLs 
avoid these drawbacks, but are susceptible to radiation damage when flux and energy are 
high (e.g. in bright-held XRM)^^. However, eq. 1 also indicates an alternative route to 
large NAs by using lenses of inferior refractive medium but with drastically miniaturized 
dimensions (i.e. small R and lenslet thickness T). 

We therefore propose Si-based objectives produced using the same nano-manufacturing 
methods as for nano-focusing ID CRLs, which are both cost-effective and capable of apex 
radii as small as 1 pm^"^. 2D focusing can be achieved by ‘interdigitating’ the ID CRLs 
such that horizontal and vertical lenslets alternate^^ as shown in Fig. 1. This reduces 
astigmatic aberration compared to sequential chips and potentially enables lens-by-lens op¬ 
timization of the NA through aberration-corrected^^, adiabatic^® and kinoform geometries^^. 
Ultimately, we show theoretically that the technical benehts of utilising Si as a refractive 
medium overcomes its inherent disadvantages in terms of refractive performance to yield 
higher performing lenses at hard energies. 

The formalism^® uses the ray-transfer matrix approach^®^^^ to describe the cumulative 
effects of the individual lenslets in the CRL and predict the aberration induced by misalign- 
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ment. At its core is the general expression for the transmission function I/Iq. In the ID 
case of axisymmetric imaging lenses, J/Jq describes the attenuation of a ray originating from 
a Gaussian source with radial position and angle Wg by a CRL comprising N identical 
parabolic lenslets of web thickness Tq and linear attenuation coefficient /i: 


I 

To 


Wg) oc exp {—NTofi) x exp 
■-(w^ -'yvgf 


X exp 


2^2 
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Eq. 2 is a product of three exponential terms: a constant absorption factor, a Gaussian 
with standard deviation describing the vignetting at the sample plane (i.e. reduction 
in image brightness towards the periphery) and another Gaussian describing the angular 
acceptance of the lens whose standard deviation (Jq is half the NA. We can therefore op¬ 
timize the NA in terms of N and R (and consequently the sample-objective distance di) 
for a given material, magnihcation and x-ray energy. This then enables the comparison of 
the best practically-achievable performance of GRLs of different materials geometries and 
conhgurations across the hard energy regime. Fig. 2 shows the optimum NA and the total 
integrated image intensity from 15 to 75 keV for a selection of imaging objectives: Axisym¬ 
metric lenses of Be and A1 (both R = 50 /rm), and interdigitated lenses of PMMA {R = 10 
pm). Si (R = 20 pm and R = 5 pm). The advantages of small-R interdigitated Si GRLs 
are evident above ~30-35 keV, where they can be more efficient than the PMMA, Be or A1 
GRLs. 

We produced a prototype interdigitated lens comprising two chips of 20 ID Si lenslets 
with 20.5±0.016 pm apex radii. The chips were produced through the standard contact 
UV-lithographic and deep reactive ion etching process^^ to a depth of 350 pm, laser-cut 
from the wafer and assembled on a steel gauge block using a micromanipulator and an 
optical microscrope. The 2D lens was tested in a full-held XRM at beamline ID06 at the 
European Synchrotron using incoherent x-rays of energy 17 keV (wavelength A = 0.7293 A). 
The microscope had a magnihcation ratio of Ad=11 and was intentionally dihraction-limited 
by the objective and oversampled by the GGD detector to reduce the contribution of the 
detector optics to the optical transfer function and resolution. 

A direct and straightforward measure of the resolution and aberration of the microscope 
can be made by inspecting the image of an absorbing Siemens star with radial features 
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FIG. 2. NA and total integrated image intensity {Itot/Io) of axisymmetric Be and A1 lenses com¬ 
pared to interdigitated polymer and Si lenses in the hard x-ray regime. The number of lenses (and 
hence focal length) is optimized for each CRT and energy to produce the best compromise of NA 
and Itot/Io- Details of the optimization routine and the CRT geometries are given in ref 18 



FIG. 3. Image of resolution chart at AI=11. The image has been post-processed by background 
subtraction and flat-held correction 

ranging from 5 fim to 50 nm (Fig. 3). The image shows a rapid degradation of contrast 
between 500-200 nm as well as horizontal striations of varying intensity. We attribute the 
contrast variation across the image to a combination of partial coherence of the incident 
beam and small shape errors in the lenslets. 
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FIG. 4. Comparison of (a) heterogeniety in the flat-held residual and (b) surface topology of a 
single lens 


Imaging the flat wavefront (i.e. without the sample) yields a magnihed intensity distri¬ 
bution related to the vignetting function in eq. 2. Fitting a 2D Gaussian to this distribution 
and inspecting its broadness and residual provides a means to compare the attenuation pro- 
hle of the CRL to theoretical expectations and to characterise shape errors in lens profile. 
For this lens, the Gaussian intensity function had standard deviations of 19.06 pm and 18.53 
pm in the horizontal and vertical directions respectively, which compare favourably with the 
theoretically-calculated values of 18.57 pm and 18.44 pm. The ~2.8% difference between the 
horizontal and vertical values suggest a slight astigmatism. The residual from this fit (Fig. 
4(a)) shows diffuse horizontal striations across the held of view. The surface topography of a 
single lenslet measured by atomic force microscopy^^ shows ridges of similar location, aspect 
ratio and frequency, suggesting that the diffuse striations and irregular image contrast may 
originate from shape errors incurred during etching (Fig. 4(b)). 

A more quantifiable measure of optical performance comes from the optical transfer 
function (OTF), 'Hifxjfv), represented in a convolution with the intensity distribution at 
the sample I/Io{xs,ys) to result in the intensity distribution at the detector I/Io{xci,ydY^ 
(where denotes the inverse Fourier transform): 
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I{xd,yd)=J^ Hxs,ys) (3) 

For a non-aberrated axisymmetric CRL, (the modulation transfer function, MTF) 
is a Gaussian function in the frequency (/r) domain: 

|W(/H)t = exp(^-^) (4) 

It then follows that by imaging a sample whose Fourier transform (FT) is known a 
priori, we can determine the MTF and calculate (Ja and thus NA. To this end, we measured 
a high-resolution image of a square grid of absorbing gold with band width W = 5 /rm 
and period P = 10 /im. The FT of such a grid is well approximated by a series of 2D 
laplacian distributions with broadness ai- Thus, including l-i and a linear background: 
B^/r) = (foo + bifn), the spatial frequency distribution of the image, G^/r) is: 


G{fR) = B{fn)+H{fR) 


N 


•s^ . ( 27rnhF \ 

X 2^ smc I ^ 1 exp 

n=-N ' 


27 rn 

P 


(XL 


(5) 


Fig. 5 shows the results from htting the horizontal and vertical FTs of the grid pattern 
with Eq. 5. This yielded values for Oa of 5.94 /rrad and 5.33 /rrad in the vertical and 
horizontal directions, which compare favourably to the theoretical predictions of 5.39 /rrad 
and 5.36 /irad. The cutoff frequency /o was solved numerically from the £t results^^ and used 
to calculate resolution values of 230 and 280 nm for the vertical and horizontal directions. 
The difference between the (Ja and /o values are also consistent with astigmatism, which 
was measured directly as the difference in focal position A/ where /o is maximized in the 
horizontal and vertical directions. The 8 mm difference can be attributed to an angular 
misalignment (i.e. non-perpendicularity of the two chips) of ~ 0.75°, which is related to the 
lenslet focal length / and misalignment angle Q according to eq. 6: 


= fx - fy ^ 


2/sin (26*) 


f 
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( 6 ) 
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While small astigmatisms are typically trivial due to the large depth of held for CRL- 
based XRMs®, we note that any predictable aberration could be remedied by the addition 
of corrective focusing elements^^. This highlights an important advantage of interdigitated 
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FIG. 5. Horizontal and vertical Fourier tranforms of the grid image and their respective fits with 
the theory in eq. 5 

geometries in the ability to incorporate optimized and aberration-corrected geometries. One 
could maximise the NA beyond that theorized in Fig. 2 by optimizing a) refractive power, 
by progressively altering the lens prohle through the CRL (i.e. adiabatic geometry^®) and b) 
transmission, by adopting kinoform or Fresnel geometries^^. We estimate that improvements 
of 1.5-2 X to the NA and 2-4 x to the transmission are within reach, constituting a signihcant 
step towards high-speed lenses for dynamic XRM. We also note that this technology could 
make even higher x-ray energies accessible through electroplated Ni lenses^®. 

In conclusion, we have shown that due to their potential for miniaturization and short 
focal lengths. Si-based interdigitated optics can in principle outperform conventional, De¬ 
based optics above 30 keV. To validate this prediction, we demonstrated a functioning 2D 
silicon objective for use in a full-held XRM at hard energies. The results are promising; 
showing acceptably low aberration and performance close to theoretical expectations. By 
harnessing the potential for miniaturization and aberration-correction of Si-based nano- 
focusing lenses, one could signihcantly improve the focusing power, transmission and NA 
at hard energies, potentially realizing hard-XRM with sub-100 nm resolution. Furthermore, 
the improved focal power available to these lenses would enable shorter imaging distances 
that are both more mechanically stable and practically achievable in modest synchrotron 










beamline hutches or laboratories. 
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